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a ■ Abstract: 

O ■ A new lattice model is presented for correlated electrons on the unrestricted 4 -dimensional 

O 

electronic Hilbert space (E>n=iC'^ (where L is the lattice length). It is a supersymmetric general- 
ization of the Hubbard model, but differs from the extended Hubbard model proposed by Essler, 
Korepin and Schoutens. The supersymmetry algebra of the new model is superalgebra gl{2\l). 
The model contains one symmetry-preserving free real parameter which is the Hubbard inter- 
action parameter U, and has its origin here in the one-parameter family of inequivalent typical 
4-dimensional irreps of gl{2\l). On a one-dimensional lattice, the model is exactly solvable by 
the Bethe ansatz. 
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The Hubbard model and the t-J model, both models of correlated electrons on a lattice, and 
exactly solvable in one dimension, have been extensively studied due to their promising role in 
theoretical condensed-matter physics and possibly in high- Tc superconductivity. The t-J model 
is a lattice model on the restricted 3'^-dimensional electronic Hilbert space (2)n=iC^ (throughout 
the paper, L is the lattice length), where the occurrence of two electrons on the same lattice site 
is forbidden. With the special choice of parameters: t = 1 and J = 2, the t-J model becomes 
supersymmetric with the symmetry algebra being the superalgebra ^^(211) ^. In ^, 
Essler, Korepin and Schoutens (EKS) proposed a model, the so-called extended Hubbard model, 
of correlated electrons on the unrestricted 4^-dimensional electronic Hilbert space (8)^=10^ . This 
EKS model, which allows doubly occupied sites and combines and extends some of the interesting 
features of the Hubbard model and the t-J model, is exactly solvable in one dimension and has 
gl{2\2) supersymmetry. 

In this Letter, we propose another direction of generalization of the Hubbard model. Specifi- 
cally, we propose a new model on the same unrestricted 4^-dimensional electronic Hilbert space 
(8>^^^C^, but with quite different interaction terms from the ones in the EKS model. Our model 
has gl{2\l) supersymmetry and contains one symmetry-preserving free real parameter which is 
exactly the Hubbard interaction parameter U; this real parameter U has its origin here in the 
one-parameter family of inequivalent typical 4-dimensional irreps of gl{2\l). The model can nat- 
urally be regarded as a modified Hubbard model with additional nearest-neighbor interactions 
and is again exactly solvable on a one dimensional lattice. The exact solvability of our model 
in one dimension comes from the fact that as an abstract dynamical model it is derived from a 
g/(2|l)-invariant rational R-matrix which satisfies the (graded) quantum Yang-Baxter equation 
(QYBE). 

It seems that only a 5/(2|l)-symmetric lattice model on the unrestricted 4^-dimensional elec- 
tronic Hilbert space could be a natural candidate for the lattice analogue of N=2 superconformal 
field theory, of which the gl{2\l) = osp{2\2) algebra defines the underlying symmetry, and which 
is a critically fixed point of the N=2 supersymmetric Landau-Ginzburg model |^]. This gives 
another motivation for our model. 

Let us begin by introducing some notation as in Q. Electrons on a lattice are described 
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by canonical Fermi operators Cj^g- and satisfying the anti-commutation relations given by 
{'^Icr' = ^ij^ar, where i,j = 1,2,---,L and cr, r =1, [. The operator Cj^o- annihilates an 
electron of spin a at site i, which implies that the Fock vacuum |0) satisfies q^o-IO >= 0. At a 
given lattice site i there are four possible electronic states: 

io), ir)i = 4io), ii), = 4jo), in). = 44^10). (1) 

By rii^cr = c\ ^Ci^cr we denote the number operator for electrons with spin a on site i, and we 
write rii = rii^^ + rii^i. The spin operators S , , S^, (in the following, the global operator O 
will be always expressed in terms of the local one Oi as O = J2i=i ™ dimension) 

form an sl{2) algebra and they commute with the hamiltonians that we consider below. 

In what follows, we only consider periodic lattice of length L. The well-known Hubbard 
model hamiltonian takes the following form: 



<ij> (7=T,i 



<ij> 



(^i,T - -\) + (^i,T - ^)K-a - \) 



(3) 



where < i,j > denote nearest neighour links on the lattice. It contains the hopping term for 
electrons and an on-site interaction term for electron pairs (coupling U). 

In m, Essler et al proposed a supersymmetric generalization of the Hubbard model. The 
supersymmetry algebra in their model is gl{2\2). We present here another supersymmetric 
generalization of the Hubbard model. The hamiltonian for our new model on a general d- 
dimensional lattice reads 

<i,j> <i,j> o-=T,i 

+(1 + + E E + 4^^Q,,)(ni,_^ + nj- _,) 

+^— E + (4) 

<i,j> 
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As will be seen, the super symmetry algebra underlying this model is gl(2\l). Remarkably, the 
model still contains the parameter U as a free parameter without breaking the super symmetry. 
Also this model is exactly solvable on the one dimensional periodic lattice, as is seen below. 
Throughout this paper, we will restrict U to the range U > —1. 

The hamiltonian (^) is obviously invariant under spin-reflection Cj^| <-> Cj^|. It can be viewed 
as an extended Hubbard model with additional nearest-neighbor interaction terms in a different 
fashion from the one proposed in . The physical nature of the additional terms is the following. 
The second term is nothing but a pair-hopping term. The third and fourth terms are bond- 
charge two-body and bond-charge-charge three-body interaction terms, respectively. And the 
last term is just a chemical potential. Clearly one can add to the above hamiltonian an arbitrary 
chemical potential (coefficient //) term and an external magnetic field (coefficient h) term 

hJ2i{''^i,l ~ '>T'i,^)j which commute with H^{U) but break its gl{2\l) supersymmetry. 

An interesting feature of our model is the discontinuity at U = . When f7 — > 0"^, the 
hamiltonian (^) reduces to 

/7Q(o+) = - E E (4,.c„. + 4,,c,.) 

<i,j> (T=T,i 

+2 E E (4,aCj,<T + 4^Ci,^)(ni,-<x-nj,_<,)2 

<i,j> (7=Ta 

+ E (ni + rij), (5) 

<i,j> 

containing a hopping term plus a bond-charge interaction term (up to a chemical potential). 
Whereas as f7 ^ 0~ , only a hopping term (and a chemical potential) survives in the hamiltonian 
®- 

Our local hamiltonian H^-{U) does not act as graded permutation of the electron states 
dl]) at sites i and j, in contrast to the hamiltonian in |^]. Nevertheless, it is supersymmetric, 
and the global hamiltonian commutes with global number operators of spin up and spin down, 
respectively. There are four supersymmetry generators for H^iU): Q|, Q|, Q^, and Qj with 
the corresponding local operators given by 

Qi,-] = -Vaui^id^^ + Va + 1 (1 - ni,|)ci,|, 

Qi^l = -^/aui^-^Ci^i + \/cyTT{l - ni_|)ci,| (6) 
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where < arg\/Z < tt, Z = a ot a + 1, and a > or a < —1 is the inverse of U: 



a 



U' 



(7) 



These generators, together with S, S"^, and two others (-El + -^1 E^, defined below), form 
the superalgebra gl{2\l). To prove this, we denote the generators of gl{2\l) by E^, /3, 7 = 1, 2, 3 
with grading [1] = [2] = 0, [3] = 1. In a typical 4-dimensional representation of gl{2\l), the 
highest weight itself of the representation depends on the free parameter a, thus giving rise to 
a one-parameter family of inequivalent irreps . Choose the following basis 



|4> 








, |3> 





1 

V"/ 



1 










(8) 



with |1), |4) even (bosonic) and |2), |3) odd (fermionic). Then in this typical 4-dimensional 
representation, E^ are 4x4 supermatrices of the form 



El = \2){3\, e! = \3){21 El = -\3){3\-\4){Al ii;^ = -|2)(2| - |4)(4|, 
El = ^|1)(2| + V^TT|3)(4|, El = V^|2)(l| +^/^TT|4)(3|, 



E' 



a|l)(3| + V^TT|2)(4|, Ef 



|3)(1|+V^TT|4)(2|, 



El = a |1)(1| + (a + 1) (|2)(2| + |3)(3|) + (a + 2) |4)(4|. 



(9) 



For a > 0, 



(10) 



and we call the representation unitary of type I. For a < —1, we have 



(11) 



and we refer to the representation as unitary of type II. In this paper, we are interested in 
these unitary representations. For a description and classification of the two types of unitary 
representations, see |8|. 
Further choosing 



|4)^|0), |3)^|T), |2) = U), ii) = in) 



(12) 



and with the help of the following identities, 

io)(oi + u)ai + iT)(Ti + in>mi = i, 

I Ti)(Ti I = n^ni, 
I TXT I = n^-n^ni, 

\ i){i\=ni-n^ni, (13) 

one can easily establish that 

ni = ia + 2)- [e'I),, ni,i = + m,^ = [e^) , + 1, (14) 

Ql^ = {eI)1 , = [eI)^, qI = [Ef)l , Q,,T = [Ef), . (15) 

The verification that the hamiltonian H^{U) commutes with all nine generators of gl{2\l) is 
just a straightforward calculation. Eq.(p!^) makes it clear that H^{U) commutes with the global 
number operators of spin up and spin down, respectively. 

The model is exactly solvable in one dimension by the Bethe ansatz. To show this, we first of 
all show that the local hamiltonian ff*°^-_(_-|^([/) on the one dimensional lattice is actually derived 
from a 5/(2|l)-invariant rational R-matrix which satisfies the (graded) QYBE. To this end, let 
C/q[5/(2|l)] be the well-known quantum (or q) deformation of gl{2\l) and V be the C/q[5/(2|l)]- 
module with highest weight (0, 0|a), which affords the g-deformed version of the one-parameter 



family of the inequivalent typical 4-dimensional irreps IC]. Without loss of generality, we 
assume q to be real. We also assume q to be generic, i.e. it is not a root of unity. For a > 
or a < —1, the module V is unitary of type I and of type II, respectively, and thus the tensor 
product V V is completely reducible. We write V V = where Vi, V2 and 

V3 are f/Jfi(/(2|l)]-modules with highest weights (0,0|2a), (0,-l|2a + 1) and (-l,-l|2a + 2), 
respectively Q, and let Pk, k = 1,2,3 be the projection operator from V onto V^. The 
trigonometric R-matrix R{x) G End{V (d) V), which satisfies the (graded) QYBE, 

(/ R{x)){R{xy) «)/)(/ ® R{y)) = {R{y) ® I){I ® R{xy)){R{x) ® I), (16) 

was given in IC, 11, in the form 



^(^) = + ^2 + ^W^3. (17) 

1 — xg^" X — 



Note, however, that q and a are both free parameters which do not enter the (graded) QYBE. 
Setting X = (f and taking the q = \ hmit, one gets the corresponding rational R-matrix (which 
also satisfies the (graded) QYBE) 



2" p(0) p(0) 



+ 2a + 2 -(0) 



(18) 



6' + 2a ^ ^ e-2a-2 
where ■, k = 1,2,3, are classical {q = 1) versions of P^, i.e. projection operator from 
y(0) yio) onto vl'^\ with ^nd 14^°^ being the g = 1 versions of V and V^, respectively. 
Note that V^^^ and V^^'^^ are actually (//(2|l)-modules, and V^^^ affords the representation (P). 
The projectors pjf^ can easily be evaluated: 



6(0) 



P. 



I -P. 



6(0) _ p(0) 



3 (19) 

where and 1^'^), = 1,2,3,4, form the symmetry adapted bases for the spaces and 
V^'^\ respectively. Note that R'{0) = I. We now compute |^'^) and {"^l), k = 1,2,3,4. By 
means of the matrix representation (0), one can show 



|1)0|1), 
1 

72 
1 

71 



(|2)0|1) + |1)®|2)), 
(|3)0|1) + |1)®|3)), 

[V^TT(|4) |1) + |1) |4)) + V^(|2) |3) - |3) |2))], 
a(|4) |1) + |1) |4)) + V^TT(-|2) |3) + |3) |2))], 



1 



v/2(2a + 1) 
1 



v/2(2a + 1) 
1 



— (|2)0|4) + |4)®|2)), 



(|3)0|4) + |4)®|3)), 



1 

72 

|4) |4) 



(20) 



which are easily seen to be orthonormal, so that 



{ml\ = {\ml))\ = A; =1,2,3,4, 

(|/?)0|^))t = (_l)[l/3>][l7>](|^))t^(|^))t^ 

(|/?))^ = (/3|, V/3 = 1,2,3,4. 



(21) 



Here stands for the grading of the state \0): = for even (bosonic) \0) and = 1 

for odd (fermionic) \[3). (Readers should keep in mind that the multiphcation rule for the tensor 
product is defined by 

(a ®b)[c®d) = (-l)WW (ac ® hd) (22) 

for any elements a, h, c and d.) 

Using the rational R-matrix (ll^) and denoting 



Rh+M = i0---i0 R'{e) 0/ • • • / (23) 

i i+1 



one may define |13| the local hamiltonian 

= =--(Pn ■ (24) 

By (0, (|20|), (^) and (^), and after tedious but straightforward manipulation, one gets, up 
to a constant, 

^S+i(f^) = -2(« + l)^M+i(«) (25) 

which implies that the local hamiltonian H^-_^_^{U) is indeed derived from the 5(/(2|l)-invariant 
rational R-matrix which satisfies the (graded) QYBE. (Note that the identity ( |25|) also indicates 
that H^{U) commutes with all the nine generators of gl{2\l) since the rational R-matrix W{6) 
is a (7/(21 1) invariant.) 

Now the exact solvability on the one-dimensional periodic lattice of our model is seen as the 
following four steps. Step 1: The hamiltonian H^{U) is self-adjoint and thus is diagonalizable. 



Step 2: Relation (25) immediately makes it clear [|13[ that on the one dimensional periodic 
lattice the global hamiltonian H^{U) commutes with the transfer matrix t{6) constructed from 
the rational R-matrix (|l8|) (see c.f. |]l^] for the standard definition of the transfer matrix), 
for any value of the parameter 6. Step 3: Using the fact, established in the rational case, 
that R^{6)^ = R^{9), where R^{6) = PR^{6) and P is the graded permutation operator of the 
electron states (||) , one may show that the transfer matrix t{0) is self-adjoint and consequently 
diagonalizable for any given parameter 6. (This result is in fact established for real but should 
also be valid for all complex 9 by using analytical continuation arguments.) We remark here 
that the results in this step are actually quite general: they are valid for any (other) rational 



R-matrices arising from unitary representations of any (other) quantum super algebras. Step 4- 
It can easily be shown that [t{6),t{9')] = 0, \f9, 9' and thus t{9) is diagonalizable simultaneously 
for all 9. Summarizing the above four steps, one sees that the hamiltonian H^{U) satisfies the 
standard requirement for a model to be exactly solvable by the Bethe ansatz. This completes 
the proof for the exact solvability in one dimension of our model. The details of solution of the 
model is deferred to a separate publication. 
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